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摘 要 本 文 基于 多 源 X 射线 层 析 成 像 技 术 及 联合 代数 重建 算法 (SART) 首次 重 构 出 黏 性 B 类 颗粒 流 化 时 局 部 气泡 的 
三 维 形态 , 并 结合 压力 脉动 分 析 技 术 获 得 了 丰富 的 全 局 流 态 化 信息 . 实验 结果 表明 ,， 黏 性 力 的 增加 强化 了 气泡 间 的 合并 作 
用 , 导致 气泡 的 尺寸 上 升 , 频率 和 速度 下 降 ， 当 床 温 超过 35°C 后 ， 节 消 状 态 被 激发 ， 并 与 自由 流 化 状态 交替 产生 . 而 且 ， 
黏 性 力 越 大 , 节 涌 在 流 化 过 程 中 所 占 的 时 间 比 例 越 长 . 
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Experimental Study on the Fluidization Dynamics of Cohesive 
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MA JiLiangl2 LIU Dao-Yinl LIANG Cail CHEN Xiao-Ping! 


(1. Key Laboratory of Energy Thermal Conversion and Control of Ministry of Education, School of Energy and Environment, 
Southeast University, Nanjing 210096 China; 
2. Key Laboratory of Coal Gasification and Energy Chemical Engineering of Ministry of Education, East China University of 
Science and Technology, Shanghai 200237, China ) 


Abstract This paper first reconstructed three dimensional version of local bubbles in a fluidized bed 
processing cohesive Geldart B particles with a combination of X-ray tomography and Simultaneous 
Algebraic Reconstruction Technique (SART). Additionally, global fluidization dynamics were also 
obtained using in-bed pressure fluctuation analysis. The experimental results show that the preserve 
of cohesive force facilitates bubble coalescence, thus increasing the bubble size and decreasing the 
bubble frequency and bubble rising velocity. When the bed temperature exceeds 35°C, slugging takes 
place alternatively with normal fluidization. Moreover, the duration of slugging sharply increases with 
the cohesive force. 
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流 化 床 反应 器 因 其 高 效 的 传 热 传 质 与 连续 处 理 ” 流 态 化 特性 研究 起 步 较 早 ， 最 初 集中 在 湿 颗 粒 的 临 
大 量 颗粒 的 能 力 ， 在 不 同 工 业 领 域 中 得 到 了 广泛 的 ，” 界 流 化 特性 中 3， 近年 来 ， 研 究 重 点 逐渐 向 流 型 转 
应 用 . 伴随 着 流 化 床 处 理 对 象 与 手段 的 多 样 化 , 颗粒 ”化 6 乳化 相 流 动 形态 四 及 流 化 失效 判定 四 等 方 
间 经 常 出 现 末 性 力 (如 干燥 、 造 粒 、 涂 层 等 过 程 ), 导 i 迁移 . 需要 指出 的 是 , 从 实验 角度 出 发 , 在 颗粒 间 
致 颗粒 的 流动 性 降低 ， 流 态 化 行为 发 生 改 变 ， 从 而 ”均匀 引入 液 桥 力 的 难度 较 大 ， 导 致 条 性 力 易 呈 各 向 
影响 流 化 床 内 的 反应 效率 以 及 反应 器 的 稳定 运行 。 “异性 ， 降 低 结 果 的 准确 性 . 针对 此 不 足 ，Shabanian 
因此 ， 揭 示 素 性 力 对 流 态 化 特性 的 影响 规律 和 内 在 。 ”等 采用 “多 聚 物 涂 层 " 法 引入 固 桥 力作 为 接触 式 黏 性 


机 理 ， 对 流 化 床 反 应 器 的 设计 和 运行 具有 重要 的 指 。 力 bo, 并 系统 地 研究 了 不 同 务 性 力作 用 下 流 化 状态 
导 意 义 ， 从 鼓 泡 流 化 向 清流 流 化 转变 特性 的 差异 552 
黏 性 力 的 种 类 众多 ， 本 文 着 重 探讨 接触 式 攻 性 由 上 可 知 ， 已 有 研究 多 集中 在 黏 性 力 对 乳化 相 


力 ， 即 在 颗粒 接触 后 产生 的 吸引 力 趾 ， 其 中 以 液 桥 ”流动 行为 的 影响 规律 上 , 对 于 气泡 相 的 研究 较 少 , 而 
力 外 和 固 桥 力 回 最 具 代表 性 . 对 于 液 桥 力 存在 下 的 ”气泡 特性 恰恰 是 影响 鼓 泡 流 化 床 运行 中 传 热 传 质 的 
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重要 参数 , 此 外 ， 流 态 化 在 高 共性 力作 用 下 的 失效 测 点 构成 . 为 了 避免 X 射线 辐射 对 实验 人 员 健 康 的 
机 理 ， 尤 其 是 失效 前 后 动力 学 参数 的 变化 规律 尚未 ”影响 , 所 有 实验 设备 均 布 置 在 由 5 mm 铝板 围 成 的 
见报 道 . 房间 中 . 实验 时 通过 电脑 和 数据 采集 卡 实现 对 系统 
本 文采 用 实验 手段 研究 黏 性 B 类 颗粒 在 三 维 鼓 。 启 停 和 实验 数据 采集 的 控制 . X 射线 强度 数据 的 采 
泡 流 化 床 中 的 流 化 特性 : 利用 “多 聚 物 涂 层 ”法 引入 集 频率 为 2500 Hz, 采集 时 长 1 min; 压力 脉动 信号 
接触 式 猪 性 力 ， 并 通过 调节 床 温 实现 对 黏 性 力 的 精确 。 的 采集 频率 为 1250 Hz, 采集 时 长 10 min. 
深 制 ; 基于 荷兰 代 尔 夫 特 理工 大 学 搭建 的 多 源 XX 射线 。 1.2 实验 材料 
层 析 成 像 系 统 重 构 床 内 的 局 部 三 维 气泡 形态 , 同时 结 本 实验 采用 “多 聚 物 涂 层 ” 法 引入 素性 力 9， 
合 压 力 脉动 分 析 技 术 提取 全 局 流 态 化 信息 ， 重 点 关注 ” 即 在 惰性 颗粒 表面 均匀 包 衰 一 层 多 聚 物 ， 利 用 多 聚 
流 化 失效 前 后 流体 动力 学 特征 的 变化 规 得 


物 茜 度 对 温度 的 高 敏感 性 在 颗粒 间 引 发 固 桥 力 即 接 

] 实验 到 钥 式 铬 性 力 。 该 方法 的 优点 是 重复 性 好 ， 医 性 力 分 
实验 系统 及 方法 布 均匀 ， 操 作 温 度 低 . 实验 采用 直径 0.6 mm, 密度 
和 2500 kg/m? 的 臻 璃 珠 作为 “基底 颗粒 "， 其 临界 流 


图 1 所 示 为 三 维 流 化 床 实验 系统 图 . 流 化 床 本 ”化 风速 Das 为 0.325 m/s， 属 于 Geldart B 类 颗粒 ; 
体内 径 140 mm， 高 1800 mm， 采 用 金属 烧结 板 布 ”采用 聚 乙 基 丙 烯 酸 酯 (PEA) 与 聚 甲 基 丙烯酸 甲 酯 
风 . 实验 通过 电 加 热 器 及 PID 控制 器 调节 流 化 风 的 PMMA) 的 混合 物 作为 多 聚 物 涂 层 . 涂 层 过程 在 滚 
温度 ， 从 而 达到 控制 床 层 温度 及 颗粒 间 黏 性 力 的 目 。 简 式 包 衣 机 内 完成 ， 操 作 步 又 已 在 文献 [10] 中 详细 
的 。 床 层 温度 由 布置 在 距 布 风 板 110 mm 处 的 热电 ”介绍 . 图 2 为 涂 层 前 后 颗粒 表面 的 SEM 图 . 本 文中 ， 
偶 进行 监测 。 在 距离 布 风 板 30 mm、110 mm、190 多 聚 物 涂 层 厚度 约 为 10 hm。 

mm 和 270 mm 处 分 别 布置 了 Kistler 7261 型 高 敏 
压力 传感器 对 床 内 表 压 进行 测量 . 同时 , 基于 X 射 
线 层 析 成 像 技 术 检 测 横 截 面 上 的 空隙 分 布 ， 其 原理 
是 当 X 射线 穿 过 颗粒 层 时 ， 射 线 强度 会 发 生 指数 衰 
减 : T= Joe-*, 式 中 jo 和 了 分 别 为 穿 过 颗粒 层 前 
后 X 射线 的 强度 ,x 为 颗粒 层 长 度 ,/ 为 衰减 系数 . 
实验 时 ,通过 测量 每 条 X 射线 强度 的 衰减 率 即 可 获 


一 、 


得 该 条 射线 经 过 的 颗粒 层 长 度 .本 实验 采用 三 个 X 图 2 颗粒 涂 层 前 后 的 表面 特性 
射线 源 互 呈 120° 角 布 置 在 流 化 床 周围 ; 每 个 射线 源 (a) 涂 层 前 ,(b) 涂 层 后 
对 应 一 组 探测 器 ， 布置 在 流 化 床 外 | 与 身 线 源 相对 Fig. 2 Surface morphology of particles before and after 


coating (a) coated (b) uncoated 


的 位 置 ; 每 组 探测 器 分 为 上 下 两 排 , 每 排 由 32 个 探 


1.3 数据 处 理 方法 


六 本 文采 用 的 压力 脉动 分 析 算法 主要 包括 功率 谱 
bea pt 、 密度 和 平均 频率 等 常用 方法 ， 因 此 不 再 歼 述 . 下 面 
人 将 对 气泡 形态 重 构 算法 进行 简要 介绍 ， 

四 Wn 气泡 重 构 过 程 分 为 三 步 : 1) 将 3 组 探测 器 共 192 


个 探测 点 检测 到 的 X 射线 强度 代入 标定 函数 ,转化 

为 每 条 X 射线 经 过 的 颗粒 层 长 度 ，P; 2) 将 重 构 区 
上 层 探测 面 域 ( 流 化 床 内 壁 的 外 切 四 边 形 ) 网 格 化 ， 并 计算 射线 
在 网 格 内 的 权重 因子 W， 本 文中 重 构 区 域 的 像素 分 
辩 率 是 80x80; 3) 求解 式 (1), 获得 每 个 像素 内 的 颗 
下 层 探测 面 粒 份额 a. 


热电 偶 


N 
pi= > Wijay (1) 
j=1 


1 三 维 鼓 泡 流 化 床 实验 系统 式 中 ， pi 是 射线 1 经 过 的 颗粒 层 长 度 ; Wij 是 身 线 7 
Fig. 1 Schematic of experimental setup 在 像素 用 内 的 权重 因子 即 射 线 长 度 ; N 是 图 像 内 的 
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像素 数 . 2 结果 与 讨论 

本 文采 用 Anderson 和 Kak 于 1984 年 提出 的 联 
合 代数 重建 算法 (SART) 求解 式 (1)。 算法 的 实现 过 图 5 是 不 同 笑 性 力作 用 下 , 距 布 风 板 70 mm 横 
程 已 在 文献 [13] 中 进行 了 详细 的 论述 . 本 文中 , 重 构 ”截面 上 的 三 维 气泡 重 构图 . 当 床 温 为 30"C 时 ,经 过 
所 用 数据 为 每 10 组 原始 数据 的 平均 值 , 意味 着 数据 ” 横 截 面 的 气泡 体积 较 小 、 数 量 众多 旦 存在 彼此 相连 
的 有 效 采样 频率 由 2500 帧 / 秒 降 为 250 帧 / 秒 . 这 昌 的 多 重 气泡 结构 ; 随 着 黏 性 力 的 上 升 ， 气 泡 频率 逐 
然 降低 了 重 构 结果 的 时 间 分 辨 率 , 但 却 有 效 地 减 小 渐 降 低 ， 气 泡 尺 寸 逐渐 增加 ， 意 味 着 气泡 间 的 合并 
了 由 射线 源 强度 波动 造成 的 误差 . 图 3 对 比 了 上 下 ”作用 增强 ; 当 床 温 达 到 45°C 时 , 床 内 出 现 了 横向 尺 
两 层 探测 面 在 同一 时 刻 的 重 构图 像 ， 图 中 白色 区 域 。 寸 与 床 内 径 相 同 的 大 气泡 , 流 化 状态 转变 为 “ 节 涌 ”。 
被 气泡 占据 . 将 不 同时 刻 的 重 构图 像 进 行 琶 加 , 即 可 ”气泡 尺寸 与 黏 性 力 间 的 正比 关系 证 明了 随 着 黏 性 力 
获得 经 过 探测 面 的 增加 ,气泡 相 与 乳化 相 之 间 的 传 质 逐 渐 亚 化。 


中 红色 对 象 为 气泡 。 
上 下 两 层 探测 面 上 出 现 的 时 间 差 估算 该 气泡 的 上 升 


的 三 维 气 泡 图 像 ， 如 图 4 所 示 ， 图 
包 . 时 


构 后 ， 可 以 基于 同一 气泡 在 


速度 . 同 理 , 也 可 以 针对 三 维 气泡 提取 尺寸 、 形 状 、 
频率 、 分 布 等 特征 参数 ， 


(b) 


图 3 上 下 两 层 X 射线 平面 重 构 结果 对 比 
(a) 上 层 探测 面 ,，(b) 下 层 探测 面 


Fig. 3 Reconstructed results for upper and lower planes 
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图 4 上 下 两 层 探测 器 3D 气泡 重 构 对 比 
(a) 上 层 探测 面 ; (b) 下 层 探测 面 
4 Reconstructed bubbles for the upper and lower planes 


(a) upper plane (b) lower plane 


1490 ys 70140 
y/mm x/mm 
图 5 距 布 风 板 70 mm 处 的 三 维 重 构 气 泡 形态 随 黏 性 力 的 变化 

规律 , 静止 床 高 : 300 mm, 流 化 风速 : 2.5Cmf 


Fig. 5 Effects of cohesive force on the bubbles reconstructed 
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at 70 mm from gas distributor. Static bed height: 300 mm, 
fluidization gas velocity: 2.5Umf 


图 6 所 示 为 不 同 黏 性 力作 用 下 各 高 度 横 截 面 上 
的 气泡 分 布 热 图 . 图 中 各 像素 的 颜色 表征 横 截 面 上 
每 一 点 被 气泡 占据 的 概率 . 在 靠近 布 风 板 的 区 域 (图 
6a, b)， 随 着 黏 性 力 的 增加 ,气泡 的 合并 增强 ,， 尺 才 
增 大 ， 因 此 床 截 面 中 心 区域 的 气泡 高 概率 区 域 逐 渐 
六 张 . 随 着 测量 高 度 的 提升 ， 气 泡 概 率 的 变化 特 怕 
与 句 性 力 紧密 相关 : 当 和 素性 力 较 低 时 (30"C)， 气 泡 
的 横向 尺寸 伴随 着 气泡 合并 逐渐 增加 ， 气 泡 高 概率 
区 的 面积 增 大 ; 当 黏 性 力 较 大 时 (45"C)， 节 涌 产 生 
的 料 栓 完 全 占据 了 横 截 面 ， 导 致 气泡 出 现 的 概率 降 
氏 ,， 因 此 测量 高 度 越 大 ， 高 概率 区 的 面积 越 小 . 
图 7 为 不 同 高 度 截面 上 的 平均 气泡 上 升 速度 随 
黏 性 力 的 变化 规律 . 随 着 攻 性 力 的 增加 ， 平 均 气 泡 
上 升 速度 逐渐 降低 . 这 与 传统 的 气泡 理论 即 自由 气 
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泡 的 上 升 速 度 与 尺寸 呈正 比 关系 [3 相 矛 盾 ， 意 味 
下 ， 床 内 存在 许多 大 尺寸 、 低 速度 
泡 ,基于 图 5 的 观测 结果 可 以 推测 ， 这 些 低速 
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气泡 是 伴随 着 节 涌 产生 的 气 栓 . 理论 上 ， 气 栓 的 上 
升 过 程 与 传统 气泡 不 同 : 由 于 气 栓 同时 受到 颗粒 重 


力 以 及 料 栓 与 边 壁 之 间 摩 擦 力 的 影响 ， 导 至 
速度 较 小 b4， 


符 随 着 屠 性 力 的 增 


其 上 升 


加 ， 节 涌 发 生 的 频 


率 提高 ， 气 栓 在 探测 到 的 气泡 中 所 占 比 例 增加 ， 导 
致 平均 气泡 上 升 速度 减 小 .需要 注意 的 是 ， 平 均 气 


泡 速度 在 床 温 超过 3 


是 床 内 节 涌 频繁 发 生 的 起 始点 ， 


5°C 后 大 幅 下 降 , 意味 着 该 温度 


和 A 30°C 40°C A 
3 
0.4 
证 | (b) 
局 
尺 0.2 
电 ) 总 等 
(a) 
0.0 
> 
黏 性 力 
图 6 不 同 高 度 上 的 气泡 分 布 热 图 随处 性 力 的 变化 规律 


(a) 70 mm; (b) 150 mm; (c) 230 mm 
静止 床 高 : 300 mm, 流 化 风速 : 2.5Dmf 


Fig. 6 Contour plots of bubble distribution as a function of 


cohesive force at different heights. Static bed height: 300 mm, 
fluidization gas velocity: 2.5Umf. (a) 70 mm; (b) 150 


mm; (c) 230 mm 
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图 7 不 同 高 度 上 的 平均 气泡 上 升 速 


床 温 /*C 


止 床 高 : 300 mm, 流 化 风速 : 2.5Umf 


度 随 竺 性 力 的 变化 规律 ， 静 


Effects of cohesive force on the averaged rising 


velocity of bubbles at different heights. Static bed height: 300 


mm, fluidization gas velocity: 2.5Umr 


度 (PSD) 随和 共性 力 的 变化 规律 . 
寺 ，PSD 强度 总 体 呈 单 峰 分 布 ， 主 
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图 8 为 不 同 高 度 处 的 压力 脉动 信号 功率 谱 密 


当 床 温 为 30°C 


主峰 外 ，PSD 
峰 , 对 应 大 气泡 的 破裂 或 多 有 
性 力 较 大 


频 在 1 Hz 左右 。 


1 线 中 还 存在 若干 能 基 较 小 的 独立 
气泡 结构 的 产生 . 当 黏 
和 ,PSD 强度 由 单 峰 分 布 转变 为 双 峰 分 布 ， 


节 涌 持续 时 间 逐 渐 延 长 ， 


图 8 


功率 谱 密度 /kPa2/Hz 


其 中 低频 峰 出 现在 0.1 Hz 左右 . 由 图 5 可 知 ， 双 峰 
分 别 对 应 节 涌 和 正常 流 化 . 而且 由 于 正常 流 化 状态 
在 持续 时 间 上 占据 主导 地 位 ， 因 此 节 涌 对 应 的 频率 
仅 有 0.1 Hz 左右 . 随 着 黏 性 力 的 增加 ， 床 内 产生 了 
大 量 的 复合 结构 气泡 [加 ， 因 此 PSD 曲线 的 频带 宽 
度 逐 渐 增 加 ; 同时 ， 曲 线 整体 开始 向 低频 方向 移动 ， 
说 明 床 内 气泡 数量 逐渐 降低 ， 其 中 ， 低 频 峰 的 峰值 


随 务 性 力 的 上 升 而 增加 ,说明 气 栓 的 尺寸 逐渐 增 大 ， 


和 | 低频 峰 


2.0 


in 


= 


So 


全 
己 


-一 30"C 
-一 407C 


”频率 /Hz 


不 同 高 度 处 压力 脉动 信号 功率 谱 密 度 (PSD) 随 柯 性 力 的 
变化 规律 , 流 化 风速 ; 2.5Uut， 静止 床 高 : 300 mm 


Fig. 8 Power spectral density of pressure signal measured at 


different heights. Fluidization gas velocity: 2.5Umt, static bed 


height: 300 mm 


图 9 为 功率 谱 密度 的 平均 频率 f。 随 苛 性 力 和 
测量 高 度 的 变化 规律 。f。 定义 为 压力 信号 功率 谱 密 


图 9 所 示 , fs 随 


压力 信号 中 处 于 低频 区 的 “气泡 通 


说 明 此 效应 


度 图 中 ， 将 信号 能 量 二 等 分 的 频率 ， 反 映 了 不 同 物 
理 现象 在 整体 信号 中 所 占 比 例 。 如 

着 共性 力 的 增加 逐渐 降低 : 在 35°C 之 前 ， 气 泡 的 
合并 作用 增强 ， 压 

道 效应 ”在 信号 中 占据 的 比例 提高 中 ， 导 致 f。 逐 
渐 降 低 ， 但 由 于 fo 的 下 降幅 度 较 小 ， 

对 压力 信号 的 能 量 分 布 影 响 有 限 ; 当 床 温 超 过 35°C 


后 , 床 内 出 现 低频 节 涌 ， 
动 ， 导致 fo 继续 下 降 . 


味 着 
上 升 


促使 能 量 分 布 向 低频 方向 移 


而 且 ， 下 降幅 度 迅 速 增加 意 


节 涌 在 流 化 过 程 中 所 占 比例 在 35°C 之 后 迅速 


201711.02565v1 


chinaXiv 


1706 工程 
2.0 T 
| 节 涌 
| —> 
N 1.6 上 
+ 
SS ! 
峙 : 
苔 | 
1.2 上 | 测量 高 度 
0 
20 和 30 35 40 45 50 
床 温 /*C 
图 9 香 性 力 及 测量 高 度 对 功率 谱 密 度 平均 频率 的 影响 , 流 化 风 


速 : 2.5Umf， 静 止 床 高 : 300 mm 
Fig. 9 Effects of cohesive force and measurement height on 
the average frequency of pressure signal. Fluidization gas 
velocity: 2.5Umf, static bed height: 300 mm 


3 结 论 
本 文 结合 X 射线 层 析 成 像 和 压力 脉动 分 析 技术 
研究 黏 性 B 类 颗粒 在 三 维 鼓 泡 流 化 床 中 的 流 化 特性 ， 
获取 了 不 同 黏 性 力 下 的 气泡 形态 演变 ， 流 化 失效 方 
式 以 及 失效 前 后 流 化 特征 参数 的 变化 规律 ， 得 到 主 
要 结论 如 下 : 
1) 当 吾 性 力 较 大 时 , 流 化 以 节 涌 的 方式 失效 , 而 
且 由 于 料 检 的 易 破碎 性 ， 床 内 交 蔡 产生 节 涌 和 正常 
流 化 ， 随 着 茜 性 力 的 增加 ， 流 化 过 程 中 节 涌 持续 时 
长 所 占 比 例 逐 渐 增 加 ; 
2) 黏 性 力 的 存在 促进 了 三 维 床 中 气泡 的 合并 ， 
导致 气泡 频率 逐渐 降低 ， 气 泡 尺 寸 逐 渐 增 加 ， 由 
于 低速 气 栓 在 气泡 中 的 占 比 与 备 性 力 呈 正比 关系 ， 
因此 床 内 的 平均 气泡 上 升 速度 随 弗 性 力 的 增加 逐渐 
降低 ; 
3) 伴随 着 节 涌 的 产生 ， 压 力 信号 的 功率 谱 密 度 
呈 双 峰 分 布 ， 低 频 峰 对 应 节 涌 ， 高 频 峰 对 应 正常 流 
化 . 节 涌 在 压力 信号 能 量 分 布 中 所 占 比例 随 黏 性 力 
的 增加 而 增加 , 导致 平均 频率 在 床 温 超过 35°C 后 迅 
速 下 降 ， 
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